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HIGHLIGHTS 


►  We  identify  two  critical  failure  modes  in  a  coated-hollow  electrode  particle. 

►  We  calculate  the  stress  fields  in  spherical  and  cylindrical  nanostructures. 

►  We  present  conditions  to  avert  fracture  and  debonding. 

►  We  discuss  the  stress  effect  on  the  lithiation  reaction. 
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In  a  novel  design  of  lithium-ion  batteries,  hollow  electrode  particles  coated  with  stiff  shells  are  used  to 
mitigate  mechanical  and  chemical  degradation.  In  particular,  silicon  anodes  of  such  core-shell  nano¬ 
structures  have  been  cycled  thousands  of  times  with  little  capacity  fading.  To  reduce  weight  and  to 
facilitate  lithium  diffusion,  the  shell  should  be  thin.  However,  to  avert  fracture  and  debonding  from  the 
core,  the  shell  must  be  sufficiently  thick.  This  tradeoff  is  considered  here  by  calculating  the  stress  fields 
resulting  from  concurrent  insertion  reaction  and  plastic  flow  for  both  spherical  and  cylindrical  hollow 
core-shell  nanostructures.  Conditions  to  avert  fracture  and  debonding  are  identified  in  terms  of  the 
radius  of  the  core,  the  thickness  of  the  shell,  and  the  state  of  charge.  The  effect  of  the  stress  on  the 
electrochemical  reaction  is  also  discussed. 

©  2012  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

High  capacity  lithium-ion  batteries  are  ubiquitous  in  portable 
electronics.  Development  of  these  batteries  is  further  motivated  by 
their  potential  use  in  electric  vehicles  [1,2],  Silicon  can  host  a  large 
amount  of  lithium,  making  it  one  of  the  most  promising  materials 
to  replace  the  conventional  carbon-based  anodes  [3],  The 
commercialization  of  silicon  anodes  is  mainly  limited  by  mechan¬ 
ical  failure  and  chemical  degradation  [4,5],  Mechanical  failure 
results  from  the  dramatic  swelling  and  stress  generated  during 
cyclic  lithiation  and  delithiation.  The  stress  may  cause  fracture  of 
the  conductive  matrix  and  active  particles,  leading  to  the  loss  of  the 
conduction  path  for  electrons,  Fig.  1(a).  As  a  result,  rapid  capacity 
fading  of  silicon  often  occurs  in  the  first  few  cycles  [6,7],  Chemical 
degradation  is  attributed  mainly  to  the  formation  of  solid 
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electrolyte  interphase  (SEI)  resulting  from  the  reactions  between 
the  electrode  and  the  electrolyte  [5],  The  formation  of  the  SEI 
consumes  active  materials  and  results  in  irreversible  capacity  loss 
in  the  first  cycle,  Fig.  1(b).  Furthermore,  the  repeated  swelling  and 
deswelling  can  cause  continual  shedding  and  re-forming  of  the  SEI 
during  cycling,  resulting  in  a  persistent  decrease  in  the  long-term 
coulombic  efficiency  [4], 

Recent  experiments  show  that  mechanical  failure  can  be  miti¬ 
gated  by  using  nanostructured  silicon  anodes,  such  as  nanowires 
[8],  thin  films  [9-12],  hollow  nanoparticles  [13],  and  nanoporous 
structures  [14],  These  structures  alleviate  the  stress  by  managing 
the  deformation  through  shape  optimization  and  geometric 
restrictions  [15],  Furthermore,  intense  efforts  are  being  dedicated 
to  stabilizing  the  SEI  layer  during  lithiation  cycles.  Current  methods 
include  coating  a  thin  artificial  SEI  layer  on  the  surface  of  silicon 
[15—21]  and  using  electrolyte  additives  to  avoid  the  continual 
consumption  of  active  materials  [22—27], 

A  recent  design  of  electrodes  involves  hollow  core-shell 
nanostructures,  which  can  enhance  the  mechanical  and  chemical 
stability  of  silicon  anodes  simultaneously,  Fig.  1(c).  Such  structures 
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a  Silicon  particle  without  a  stiff  shell 


b  Silicon  particle  without  a  stiff  shell 


C  Hollow  silicon  particle  with  a  stiff  shell 


SEI  Si 


Fig.  1.  (a).  For  a  silicon  particle  without  a  stiff  shell,  the  stress  induced  by  lithiation  and  delithiation  may  cause  fracture  of  the  active  particle,  or  debonding  between  the  active 
particle  and  the  matrix,  leading  to  the  loss  of  the  conduction  path  for  electrons,  (b).  Also  for  a  silicon  particle  without  a  stiff  shell,  the  deformation  associated  with  lithiation  and 
delithiation  may  cause  the  shedding  and  re-forming  of  the  solid-electrolyte  interphase  (SEI),  consuming  active  materials,  (c).  For  a  hollow  silicon  particle  with  a  stiff  shell,  the 
deformation  of  silicon  is  accommodated  by  inward  swelling,  so  that  electric  contact  is  maintained,  and  the  shedding  of  SEI  avoided. 


enable  silicon  electrodes  to  sustain  over  six  thousand  cycles  with 
high  coulombic  efficiency  [28],  In  a  hollow  core-shell  nano¬ 
structure,  the  shell  separates  the  electrolyte  and  the  electrode 
(core),  and  forces  the  core  to  swell  into  the  hollow  space  in  the 
interior  of  the  structure.  Because  the  shell  restricts  outward 
swelling  of  the  core,  the  continual  SEI  shedding  and  re-forming  is 
suppressed,  and  the  electric  contact  between  the  electrode  and  the 
matrix  is  maintained.  Similar  designs  based  on  encapsulating 
silicon  nanoparticles  within  hollow  carbon  shells  show  promise  for 
high-capacity  electrodes  with  long  cycle  life  of  hundreds  of  cycles 
[29-31], 

Constrained  lithiation  and  delithiation,  however,  induce  stress 
in  the  hollow  core— shell  structures.  During  lithiation,  Fig.  2(a),  the 
core  is  under  compressive  stresses,  and  the  shell  is  under  tensile 
hoop  stress.  Such  tensile  stress  can  initiate  fracture  in  the  shell. 
During  delithiation,  Fig.  2(b),  the  radial  stress  can  be  tensile  in  both 
the  core  and  the  shell.  Such  stress  may  cause  interfacial  debonding. 
To  facilitate  lithium  diffusion  and  reduce  overall  weight,  the  ideal 
stiff  shell  should  be  thin.  However,  an  extremely  thin  shell  is  unable 
to  constrain  the  outward  expansion  of  the  core.  The  insertion 


reaction  causes  large  deformation  in  the  core,  but  the  magnitude  of 
the  stress  is  kept  manageable  by  plastic  flow.  Here  we  calculate  the 
stress  field  resulting  from  the  concurrent  insertion  reaction  and 
plastic  flow  in  hollow  spherical  nanoparticles  and  nanowires  of 
silicon  coated  with  stiff  shells.  We  identify  conditions  to  avoid 
fracture  and  debonding  in  terms  of  the  radius  of  the  core,  the 
thickness  of  the  shell,  and  the  state  of  charge.  We  further  identify 
the  driving  forces  for  the  lithiation  reaction  and  discuss  the  effect  of 
the  stress  field  on  the  electrochemical  reaction. 

2.  Stresses  in  coated  hollow  spherical  particles  of  silicon 

In  a  hollow  core-shell  structure,  the  stiff  shell  forces  the  core  to 
swell  into  the  inner  hollow  space  (Fig.  3).  Before  absorbing  any 
lithium,  the  inner  radius  of  the  hollow  particle  is  A,  and  the  outer 
radius  is  B.  The  particle  is  coated  by  a  stiff  shell  of  outer  radius  C. 
Note  that  a  typical  value  of  lithium  diffusivity  in  amorphous  silicon 
is  ~10-16  m2  s-1  [32,33],  the  diffusion  time  for  lithium  through 
a  100  nm  silicon  film  can  be  estimated  as  on  the  order  of  100  s 
(t  ~  L2/D).  Such  time  is  much  shorter  than  the  time  of  a  lithiation 
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Fig.  2.  Two  potential  modes  of  failure  in  a  hollow  silicon  particle  coated  with  a  stiff  shell,  (a).  The  lithiation  of  the  silicon  particle  induces  tensile  hoop  stress  in  the  shell,  which  may 
cause  the  shell  to  fracture,  (b)  The  delithiation  of  the  silicon  particle  induces  radial  tensile  stress,  which  may  cause  debonding  between  the  core  and  the  shell. 


cycle  in  typical  electrochemical  measurements.  Thus  we  consider 
that  diffusion  in  the  core  is  fast  to  allow  for  a  uniform  distribution 
of  lithium.  To  completely  fill  the  hollow  space  upon  full  lithiation, 
the  structure  should  satisfy  the  geometric  relation: 


where  the  volumetric  swelling  ratio  0  =  t/f/Vj,  with  Vf  representing 
the  fully  lithiated  volume  and  Vi  the  initial  volume  of  the  lithium- 
free  state.  Lithiation  of  silicon  causes  a  volumetric  swelling  /?  =  4, 
so  that  A/B  =  (3/4)1/3.  [34] 

We  take  the  stress-free  and  lithium-free  state  as  the  reference 
state,  Fig.  3(a).  At  time  t,  the  hollow  particle  absorbs  some  lithium, 
retaining  its  spherical  symmetry.  We  assume  that  the  outward 
expansion  of  the  core  is  completely  constrained  by  the  stiff  shell. 
Note  that  the  elastic  deformation  of  the  shell,  typically  less  than  1%, 
is  much  less  than  the  lithiation-induced  strain  in  the  silicon  core. 
With  this  assumption,  lithiation  changes  only  the  inner  radius  a, 
and  we  neglect  any  deformation  of  the  shell,  Fig.  3(b).  We  represent 
a  material  element  in  the  reference  state  by  its  distance  R  from  the 


center.  At  time  t,  it  moves  to  a  place  a  distance  r  from  the  center.  The 
function  r(R,  t)  specifies  the  kinematics  of  the  deformation.  Each 
material  element  in  the  hollow  core  is  subject  to  a  state  of  triaxial 
stresses,  ( <rr ,  oe,  <r$),  where  a,  is  the  radial  stress  and  ag  the  hoop 
stress.  The  stresses  are  inhomogeneous,  represented  by  functions 
<7r(r,  t)  and  ag(r,  t).  The  balance  of  forces  acting  on  a  material 
element  requires  that 

9 ar(r,t)-ae(r,t)  =Q_ 

dr  r 

We  adopt  the  commonly  used  idealization  that  plastic  defor¬ 
mation  is  unaffected  when  a  hydrostatic  stress  is  superimposed. 
Furthermore,  the  yield  strength  of  the  lithiated  silicon,  oy,  is  taken 
to  be  constant  and  independent  of  the  amount  of  deformation  and 
the  concentration  of  lithium  [35].  Because  the  elastic  strain  is 
negligible  compared  to  the  lithiation-induced  strain,  the  elasticity 
of  the  hollow  core  is  neglected.  During  lithiation,  a  material 
element  of  silicon  is  under  compression  in  the  hoop  direction  and  is 
in  the  state  <rr  -  og  =  cry.  Setting  ar  -  ag  —  ay  in  Eq.  (2) 
and  integrating  over  r  with  the  traction-free  boundary  condition, 
<rr(a,  £j  =  0,  the  radial  stress  in  the  hollow  core  is  given  by 


Fig.  3.  (a).  In  the  reference  state,  a  hollow  particle  of  an  electrode  is  stress-free  and  lithium-free,  (b)  In  I 
is  accommodated  by  the  inner  hollow  space.  Outward  deformation  is  restricted  by  the  shell. 


:  particle  is  partially  lithiated.  The  deformatic 
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or  =  -2ffYlog^.  a  <  r  <  B.  (3) 

The  hoop  stress  in  the  hollow  core  is  determined  by  oy  -  ag  =  cy, 
giving 


ag  =  -<ty(21og^+  l),  a  <  r  <  B  (4) 

The  stress  fields  in  the  elastic  shell  take  the  familiar  solutions  of 
Lame  problems  [36],  with  the  radial  stress 

*=-2^(log§)(L^)(£-l).  B<r<C,  (5) 

and  the  hoop  stress 

a9  =  2<rY(log®)(c:3^3)(g+l),  B  <r  <C  (6) 

During  delithiation,  a  material  element  of  the  hollow  core  is 
under  tension  in  the  hoop  direction,  and  is  in  the  state 
Gy  -  erg  =  — ay.  Setting  crr  -  ag  =  —cry  in  Eq.  (2)  we  obtain  the  radial 
stress  in  the  hollow  core  by  integrating  over  r, 

<rr  =  2oylog-,  a  <  r  <  B,  (7) 

and  the  hoop  stress, 

=  ffy(21og~  +  l),  a<r<B  (8) 

A  comparison  of  the  stress  fields  indicates  a  jump  in  both  the 
radial  and  hoop  stresses  from  lithiation  to  delithiation,  as  shown  in 
Fig.  (4).  These  jumps  are  caused  by  neglecting  the  elastic  strain. 
Should  we  include  the  elastic  strain,  the  stresses  would  make  this 
transition  from  compressive  to  tensile  yielding  after  the  removal  of 
a  very  small  amount  of  lithium. 

Mechanical  failure  is  mainly  due  to  the  stress  at  the  interface 
between  the  core  and  shell.  Fig.  4  plots  the  evolution  of  the  radial 
stress  at  the  interface  of  the  hollow  core  during  lithiation  and 
delithiation.  The  stress  is  plotted  as  a  function  of  the  state  of  charge 
(SOC)  —  0  represents  the  lithium-free  state,  and  1  the  fully  lithiated 
state.  Given  the  geometric  condition  in  Eq.  (1 ),  the  state  of  charge  is 
calculated  by  using  the  inner  radius  a  at  a  given  state, 


State  of  charge 

Fig.  4.  Evolution  of  the  radial  stress  at  the  interface  between  the  particle  and  the  shell 
during  lithiation  and  delithiation. 


As  expected,  the  hollow  core  is  under  compression  in  the  radial 
direction  during  lithiation.  The  compressive  stress  increases  loga¬ 
rithmically  with  the  state  of  lithiation,  Eq.  (3).  Because  of  the 
triaxial  state  of  stress,  the  magnitudes  of  the  components  of  the 
stress  readily  exceed  the  yield  strength  of  lithiated  silicon  [37],  Such 
big  compressive  stresses  in  the  core  can  cause  a  large  tensile  hoop 
stress  in  the  outer  shell,  which  may  lead  to  fracture  of  the  shell. 
During  delithiation,  the  core-shell  interface  is  under  tensile  stress 
in  both  the  radial  and  hoop  directions.  The  magnitudes  of  the  stress 
components  are  maximized  at  the  beginning  of  the  delithiation, 
and  gradually  decrease  with  the  state  of  delithiation.  The  tensile 
stress  in  the  radial  direction  may  initiate  interfacial  debonding, 
resulting  in  a  loss  of  electrical  contact  between  the  active  core  and 
the  surface  coating. 


3.  Fracture  and  debonding  analysis 


We  now  analyze  fracture  of  the  coating  using  an  approach 
similar  to  that  described  in  several  recent  papers  [38—40],  We  focus 
on  fracture  of  the  coating  caused  by  the  tensile  hoop  stress  during 
lithiation.  Linear  elastic  fracture  mechanics  are  adopted.  The 
reduction  in  the  elastic  energy  associated  with  a  crack  advancing 
a  unit  area  defines  the  energy  release  rate  Gf.  Dimensional  analysis 
dictates  that  the  energy  release  rate  should  take  the  form 

Gf=zg(C-B).  (10) 

where  Es  is  Young’s  modulus  of  the  shell,  Z  a  dimensionless 
parameter  to  be  determined  by  solving  the  elastic  boundary-value 
problem,  og  the  stress  in  the  hoop  direction  given  by  Eq.  (6),  and 
C  -  B  the  thickness  of  the  surface  coating.  In  the  case  of  a  channel 
crack  in  the  thin  shell,  Z  =  2  is  a  typical  value.  Inserting  this  value  of 
Z  and  Eq.  (6)  into  Eq.  (10)  gives  an  analytical  solution  of  the  energy 
release  rate 


y^) 


C3  V 
f2B3j 


where  we  have  approximated  r  ~  B  since  B  »  C  -  B.  In  this  limit, 
Eq.  (11)  can  be  further  approximated  as 


4K) 


2  B2 
C^B' 


(12) 


The  energy  release  rate  is  proportional  to  the  square  of  the 
yield  strength  of  lithiated  silicon,  the  compliance  of  the  surface 
coating,  the  ratio  of  the  sizes  of  core  and  shell,  and  the  size  of  the 
core.  Therefore,  to  minimize  the  energy  release  rate,  it  is  desired 
to  have  a  small  and  soft  hollow  core  and  a  thick  and  stiff  surface 
shell.  However,  the  shell  also  needs  to  be  thin  to  minimize 
weight  and  to  allow  for  fast  diffusion  of  lithium  through  its 
thickness. 

We  take  the  same  approach  for  the  analysis  of  interfacial 
debonding.  Debonding  of  thin-film  electrodes  from  a  substrate 
associated  with  lithiation  cycles  has  been  observed  in  experi¬ 
ments  [41,42]  and  analyzed  in  theories  [10,43],  Here  we  consider 
that  the  driving  force  of  debonding,  the  energy  release  rate,  is 
mainly  attributed  to  mode  I  fracture  caused  by  the  tensile  radial 
stress  at  the  interface  of  the  silicon  core  and  the  stiff  shell.  We 
assume  a  pre-existing  crack  of  a  size  comparable  to  the  thickness 
of  the  shell,  C  -  B.  The  energy  release  rate  may  be  reduced  to 
a  simple  form, 
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Cd=n^(C-B).  (13) 

The  value  of  the  numerical  pre-factor  in  (13)  is  approximate, 
which  changes  somewhat  with  the  elastic  mismatch  between  the 
shell  and  the  core,  and  with  the  relative  length  of  the  crack  to  the 
thickness  of  the  shell.  Accurate  values  of  the  factor  can  be  calcu¬ 
lated  by  solving  the  boundary-value  problem  using  the  finite- 
element  method  [38],  In  this  work,  however,  a  particular  value  is 
used  to  illustrate  these  main  ideas.  The  radial  stress  during  deli- 
thiation  is  given  by  Eq.  (7),  and  thus 

0.2  /  d\  2 

Cd  =  47t^(tog-J  (C-fl).  (14) 

Here  Ee  is  the  effective  modulus  that  accounts  for  the  effects  of 
inhomogeneous  properties  of  the  core  and  shell.  A  typical  relation 
is  represented  by 


State  of  charge 


where  Ec  is  the  modulus  of  the  core,  and  Es  the  modulus  of  the  shell. 
[44]  Comparing  Eqs.  (12)  and  (14),  we  can  see  some  interesting 
differences.  Namely,  the  energy  release  rate  for  interfacial 
debonding  is  seen  to  scale  directly  with  the  shell  thickness  while 
the  energy  release  rate  for  fracture  of  the  shell  scales  inversely  with 
the  shell  thickness.  Furthermore,  the  energy  release  rate  for 
debonding  scales  much  more  slowly  with  the  core  radius  B,  than 
does  the  energy  release  rate  for  fracture  of  the  shell. 

Let  Tf  be  the  fracture  energy  of  the  shell,  and  rd  the  interfacial 
energy  between  the  core  and  the  shell.  No  preexisting  flaws  will 
advance  if  the  maximum  energy  release  rate  is  less  than  the  frac¬ 
ture  energy.  Therefore,  Gf  =  If  defines  the  critical  condition  for 
initiation  of  fracture  of  the  shell,  and  Gd  =  rd  defines  the  critical 
condition  for  initiation  of  interfacial  debonding. 

Recent  experiments  indicate  that  a  thin  layer  of  A1203  deposited 
on  silicon  anodes  by  atomic  layer  deposition  provides  excellent 
improvement  of  the  cycling  stability  of  silicon  anodes  [17,21],  Here 
we  take  the  example  of  a  hollow  spherical  particle  of  silicon  coated 
with  a  thin  shell  of  AI2O3  to  illustrate  the  above  analysis.  Repre¬ 
sentative  values  of  lithiated  silicon  and  AI2O3  are  taken,  oy  -  1  GPa 
[35]  and  Ec  =  12  GPa  for  lithiated  silicon  [45,46],  Es  =  300  GPa, 
rf  =  40  J  m-2  [47],  and  Td  =  1  J  itt2. 

We  first  construct  a  phase  diagram  to  show  the  effects  of  the 
thickness  of  the  shell  and  the  state  of  charge,  Fig.  5(a).  We  take 
a  typical  value  of  the  silicon  core  size,  B  =  200  nm  and 
A/B  =  (/3/4.  The  thickness  of  AI2O3,  C  -  B,  is  varied.  The  state  of 
charge  is  calculated  from  Eq.  (9).  The  solid  black  line  represents  the 
critical  conditions  for  fracture  of  the  shell,  corresponding  to  the 
condition  Gf  =  If,  and  the  red  dashed  line  (in  the  web  version) 
represents  the  critical  conditions  of  interfacial  debonding,  corre¬ 
sponding  to  the  condition  Gd  =  Td.  These  two  lines  delineate  an 
upper  boundary  to  the  safe  regime  in  which  no  fracture  of  the 
surface  coating  or  interfacial  debonding  would  occur.  Fig.  5(b) 
further  plots  a  phase  diagram  to  show  the  effects  of  core  size  and 
state  of  charge.  In  this  case,  the  thickness  of  the  AI2O3  shell  is  fixed, 
C  —  B  =  5  nm,  and  the  core  size  B  is  varied  while  maintaining  the 
geometric  relation  A/B  =  {/3/4.  Given  the  analytical  solutions  in 
Eqs.  (11)  and  (14),  phase  diagrams  for  other  experimental  condi¬ 
tions  can  be  readily  constructed. 

4.  Analysis  of  coated  hollow  silicon  nanowires 

One  dimensional  nanowires  and  nanotubes  are  emerging 
designs  that  can  mitigate  the  mechanical  failure  of  silicon 


State  of  charge 

Fig.  5.  Conditions  of  fracture  and  debonding  for  a  spherical  hollow  particle  plotted  in 
the  plane  of  (a)  the  thickness  of  the  shell  and  the  state  of  charge,  and  (b)  the  radius  of 
the  particle  and  the  state  of  charge. 


electrodes.  In  particular,  Y.  Cui  and  co-workers  illustrate  that 
coated  hollow  nanowires  enable  silicon  anodes  to  sustain  over  six 
thousands  cycles  with  high  coulombic  efficiencies  [28],  Such 
excellent  performance  would  ultimately  meet  the  requirement  of 
commercializing  silicon  into  electric  vehicles.  In  this  section,  we 
analyze  fracture  and  debonding  of  a  hollow  nanowire.  The  cross- 
section  is  as  shown  in  Fig.  3.  For  complete  utilization  of  the 
available  hollow  volume  upon  full  lithiation,  we  must  satisfy  the 
geometric  relation: 


where  for  silicon  0  =  4  and  A/B  =  (3/4)1/2. 

Similar  to  the  case  of  spherical  particles,  lithiation  of  silicon 
causes  a  tensile  hoop  stress  in  the  shell,  which  may  fracture  the 
coating  material.  Furthermore,  a  large  tensile  radial  stress  is 
developed  at  the  beginning  of  delithiation,  which  may  initiate 
interfacial  debonding.  We  derive  the  stress  field  in  the  Appendix. 
Since  the  deformation  of  the  core  along  the  axial  direction  is  con¬ 
strained  by  the  shell,  lithiated  silicon  is  assumed  to  deform  under 
plane-strain  conditions.  The  analytical  solution  of  fracture  energy 
release  rate  of  the  shell  during  lithiation  is  given  by 
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State  of  charge 


Fig.  6.  Conditions  of  fracture  and  debonding  for  a  hollow  nanowire  plotted  in  the 
plane  of  (a)  the  thickness  of  the  shell  and  the  state  of  charge,  and  (b)  the  radius  of  the 
particle  and  the  state  of  charge. 


^S('-I)2(c^)2('4)V-».  tn> 


In  the  case  of  C  ~  B,  Eq.  (17)  can  be  further  approximated  as 


And  the  interfacial  fracture  energy  during  delithiation  is  given 


by 

Cd 


4Tt  (Ty 


(19) 


Comparing  the  fracture  energies  in  the  spherical  particle  and  the 
nanowire,  we  find  that  the  energy  release  rate  for  fracture  of  the 
shell  is  a  factor  of  4/3  larger  for  the  nanowire  than  for  the  sphere.  In 
contrast,  energy  release  rate  for  interfacial  debonding  is  3  times 
larger  in  the  sphere  than  in  the  nanowire. 

We  now  consider  the  practical  example  of  a  hollow  silicon 
nanowire  coated  with  a  thin  layer  of  AI2O3.  A  phase  diagram 
showing  the  effects  of  the  thickness  of  the  shell  and  the  state  of 
charge  is  given  in  Fig.  6(a).  In  the  nanowires,  the  state  of  charge  is 
calculated  by 


SOC 


A2 -a2 
A2 


(20) 


The  silicon  core  size  is  taken  as  B  =  200  nm  and  A/B  —  (3/4)^2. 
The  thickness  of  AI2O3,  C  —  B,  is  varied.  The  solid  black  line 
corresponds  to  the  critical  conditions  of  fracture  of  the  shell,  given 
by  Gf  =  ff,  and  the  red  dashed  line  (in  the  web  version)  corre¬ 
sponds  to  the  critical  conditions  for  interfacial  debonding,  given  by 
Gd  —  Td-  Fig.  6(b)  shows  a  phase  diagram  demonstrating  the  effects 
of  the  core  size  and  state  of  charge.  In  this  case,  the  thickness  of  the 
AI2O3  shell  is  taken  as  C  -  B  =  5  nm  and  the  core  size  B  is  varied 
while  maintaining  the  geometric  relation  A/B  =  s/3/4.  More 
accurate  measurements  of  Tf  and  f  j  are  needed  to  produce  more 
accurate  plots.  In  the  case  where  the  fracture  energies  Tf  and  Td  are 
less  than  the  values  we  have  used  here,  the  black  and  red  dashed 
lines  in  Figs.  (5)  and  (6)  would  shift  to  lower  limits  of  the  state  of 
charge. 


5.  Stress  effect  on  the  lithiation  reaction 

We  now  identify  the  driving  forces  for  the  lithiation  reaction  of 
silicon,  namely,  the  change  of  free  energy  associated  with  the 
reaction.  Let  AGr  be  the  free  energy  of  the  lithiation  reaction  of 
silicon  when  both  the  stress  and  the  applied  voltage  vanish. 
Assume  the  silicon  anode  is  connected  to  a  cathode  by  a  conduct¬ 
ing  wire  through  a  voltage  source.  Associated  with  converting  one 
lithium  atom  and  silicon  into  lithiated  silicon,  one  electron  passes 
through  the  external  wire,  so  that  the  voltage  source  does  work 
e<J>,  where  4>  is  the  voltage,  and  e  the  elementary  charge  (a  positive 
quantity).  The  driving  force  is  further  modified  when  the  stress 
effect  is  included.  During  lithiation  of  the  coated  hollow  silicon 
particle,  lithium  atoms  migrate  into  silicon  by  overcoming  an 
energy  barrier  induced  by  the  compressive  stress  in  the  core.  Thus, 
when  one  lithium  atom  is  incorporated  into  lithiated  silicon,  the 
stress  does  work  -Qom,  where  <7m  is  the  mean  stress  in  the  core, 
and  Q  is  the  change  in  volume  caused  by  insertion  of  one  lithium 
atom.  Here  we  assume  an  isotropic  swelling  of  silicon  caused  by 
lithium  insertion  such  that  only  the  contribution  from  the  mean 
stress  is  included.  We  also  consider  small  elastic  deformation  of 
lithiated  silicon  such  that  stress  terms  of  higher  orders  can  be 
neglected  [48,49],  Combining  the  above  contributions,  we  find 
that,  when  the  reaction  advances,  the  net  change  in  the  free 
energy  is 


AG  =  AGr  -e&-  Qom.  (21) 

We  have  neglected  the  dissipation  at  the  electrolyte/electrode 
interfaces,  as  well  as  inside  the  electrodes  and  electrolytes.  In  our 
sign  convention,  a  negative  AG  drives  lithiation,  and  a  more  nega¬ 
tive  value  represents  a  larger  driving  force.  The  free  energy  of 
reaction  AGr  takes  a  negative  value.  The  sign  of  e<f>  depends  on  the 
polarity.  We  consider  <f>  to  be  positive  in  the  direction  that  drives 
lithiation.  As  expected,  a  compressive  mean  stress  in  silicon  retards 
lithiation. 

For  both  the  spherical  silicon  particles  and  nanowires,  the 
analytical  solutions  of  the  stress  fields  are  given  in  the  prior 
sections.  Inserting  the  stresses  into  Eq.  (21 ),  we  can  quantify  the 
effect  of  stress  on  the  lithiation  reaction.  We  note  that  the 
compressive  stress  is  maximized  at  the  core-shell  interface.  Thus, 
the  retarding  effect  of  stress  on  the  lithiation  reaction  is  most 
significant  at  the  interface.  The  contribution  due  to  the  stresses  at 
the  interface  is  plotted  in  Fig.  7.  In  making  this  plot,  we  have 
adopted  the  value  Q  —  1.36  x  10-29  m3  [50].  As  expected,  the 
magnitude  of  the  stress  contribution  increases  with  the  state  of 
charge.  Recall  that  the  free  energy  of  formation  of  lithiated  silicon  is 
small;  for  example,  AGr  =  -0.18  eV  for  amorphous  Li2.iSi  [51]. 
Consequently,  the  reaction  can  readily  generate  a  large  enough 
stress  to  counteract  the  electrochemical  driving  force,  stalling  the 


:  al.  /  Journal  of  Power  Sources  218  (2012)  6-14 


State  of  charge 

Fig.  7.  The  effect  of  the  stress  on  the  driving  force  for  the  lithiation  reaction  is  plotted 
as  a  function  of  the  state  of  charge. 


surface  reaction  prior  to  complete  lithiation.  In  fact,  a  lower  lith¬ 
iation  capacity  due  to  an  incomplete  surface  reaction  has  been 
observed  in  recent  experiments  [15,52,53]. 

6.  Conclusions 

We  present  an  analysis  of  fracture  and  debonding  failure  in 
coated  hollow  spherical  particles  and  nanowires  of  silicon  anodes. 
A  phase  diagram  describing  the  critical  structural  parameters  and 
operating  conditions  is  constructed,  outlining  how  to  avoid  fracture 
of  the  coating  and  debonding  between  the  coating  and  the  active 
material.  We  further  explore  the  effect  of  the  stress  field  on  the 
lithiation  reaction  of  silicon.  The  large  compressive  stress  at  the 
core-shell  interface  may  counteract  the  electrochemical  driving 
force,  leading  to  a  low  lithiation  capacity.  The  present  results  and 
associated  analysis  provide  insightful  guidelines  for  a  viable  design 
of  coated  hollow  nanostructures. 
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Appendix 

We  derive  the  stress  field  in  a  hollow  silicon  nanowire  coated 
with  a  stiff  shell.  The  geometry  of  the  cross  section  is  represented  in 
Fig.  3.  We  represent  a  material  element  in  the  reference  configu¬ 
ration,  Fig.  3(a),  by  its  distance  R  from  the  center.  At  time  t,  it  moves 
to  a  place  at  a  distance  r  from  the  center,  Fig.  3(b).  The  function  r(R, 
t)  specifies  the  deformation  kinematics  of  the  silicon  nanowire.  Due 
to  the  constraint  imposed  by  the  shell,  lithiated  silicon  is  assumed 
to  deform  under  the  plane-strain  conditions.  To  focus  on  the  main 
ideas,  we  neglect  the  elasticity  of  both  the  core  and  the  shell;  we 
model  the  lithiated  silicon  as  a  rigid-plastic  material.  Consequently, 
the  expansion  of  lithiated  silicon  is  entirely  due  to  lithiation. 
Consider  an  annulus  of  the  lithiated  silicon  between  the  radii  B  and 
r.  This  annulus  is  lithiated  from  the  annulus  of  pristine  silicon 
between  the  radii  B  and  R.  We  assume  that  the  rate  of  lithiation  is 
slow  relative  to  diffusion  of  lithium  in  the  core,  so  that  lithium 


atoms  have  enough  time  to  homogenize  in  the  core.  That  is,  the 
ratio  of  the  volume  of  the  lithiated  silicon  over  the  volume  of 
pristine  silicon  /?,  is  taken  to  be  homogeneous  and  evolves  in  time. 
Upon  full  lithiation,  0  reaches  =  4  for  silicon.  Thus, 


B2  -r2  =  -R2),  (A.1) 

This  equation  gives  the  function  r(R,  t )  once  the  functional)  is 
given.  That  is,  , 6(t)  fully  specifies  the  kinematics  of  the  silicon 
nanowire, 

r  =  \Jb2  -  (3(B2  -  R2) .  (A.2) 

The  stretches  can  be  calculated  as 


*r  =  ^  =  ^,  ^  =  n,  Az  =  l, 


We  decompose  the  stretches  by  writing 


(A.3) 


Ar  =  Xfp1/3,  Xe  =  A^V3,  Xz  =  Xl  f/3  (A.4) 

where  /?  represents  the  volume  change  due  to  the  insertion  of 
lithium,  and  the  plastic  stretch  Ap  represents  the  shape  change 
during  lithiation.  The  volume  change  has  been  assumed  to  be 
isotropic  for  amorphous  silicon.  We  can  calculate  the  strain 
components  from  the  stretches, 


eT  =  logAr,  eg  =  logXff,  ez  =  logAz  (A.5) 

The  incremental  plastic  deformation  is  given  by 


$4  =  <51ogAz  -  ^Slog/j  = 

The  equivalent  plastic  strain  is 


<54, 


/l  IB4 
V9  +  3  r4' 


We  adopt  the  flow  rule 


(A.6) 


(A.7) 


Si,  = 

1J  3  54/ V 


(A.8) 


where  sy  is  the  deviatoric  stress,  defined  as  sy  =  cry  -  (l/3)cry<5y,  and 
(Ty  the  yield  strength  of  lithiated  silicon.  Therefore, 


(A.9) 
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Or  -  Ofl  =  Sr  -  Sg  % 


V^4/B4  + 


Thus,  using  Eq.  (14),  the  energy  release  rate  for  interfacial 
debonding  takes  the  solution 


Consider  the  force  balance  of  a  material  element  in  lithiated 
silicon 


8°f  |  °V  ~  09  „ 

3r  +  r  ^  ’ 


(A.11)  References 


the  radial  stress  can  be  obtained  by  integrating  Eq.  (A.11),  it  gives 

Or  =  — ^ffylog  Vr4  +  3B^-V3B2  +  D  a<r  <B  (A.12) 

where  D  is  the  integration  constant.  With  the  traction-free 
boundary  condition,  (jr(a,  t)  =  0, 

a/3  a  ,  rjr4  +  3B4  -  s/3B2\ 

or  =  21og-  +  log  v  = - — —  ,  a  <  r  <  B 

3  l  r  a/o4  +  3B4-a/3B2/ 


The  stresses  along  hoop  and  axial  directions  are  obtained  from 
Eqs.  (A.9)  and  (A.10), 


2<jy 


(A.14) 

r  =  „  ffY _ 

Z  r  V^/B4  +  3  V3B4/r4  +  r 

Therefore  the  radial  stress  at  the  interface  r  =  B  is  given  by, 


a/3  (nl  a  ,  (2-V3)B2 

(7r  =  o^yI  21og-  +  log— z - 


(A.15) 


In  the  fully  lithiated  state  of  a  hollow  silicon  nanowire,  a  «  B. 
The  radial  stress  at  the  interface  can  then  be  approximated  by 
a  Taylor  expansion,  which  gives 


The  stress  field  in  the  elastic  shell  can  be  solved  using  the 
familiar  solution  of  Lame  problem  [36],  it  gives, 

2 a/3  /.  a\f  2 vB2  \ 

^  =  -  3  ffv  (  °®  B/  \C2  -  B2) ' 


Similar  to  the  analysis  of  a  sphere,  the  energy  release  rate  in  the 
coating  shell  takes  the  following  solution 

«=§i«)2(  c^)2('4)V*>.  *>■> 

During  delithiation,  the  tensile  radial  stress  in  silicon  at  the 
interface  becomes 
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